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Abstract. 

We present statistical analysis of acoustic emission (AE) data from tensile 
experiments on paper sheets, loading mode I, with samples broken under strain control. 
The results are based on 100 experiments on unnotched samples and 70 samples with 
a long initial edge notch. First, AE energy release and AE event rates are considered 
for both cases, to test for the presence of "critical points" in fracture. For AE energy, 
no clear signatures are found, whereas the main finding is that the event rate diverges 
when a sample-dependent "critical time" of the maximum event rate is approached. 
This takes place after the maximum stress is reached. The results are compared with 
statistical fracture models of heterogenous materials. We also discuss the dependence 
of the AE energy and event interval distributions on average event rates. 
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1. Introduction 

Fracture of heterogenous material exhibits scaling properties familiar from statistical 
physics. Scaling laws of energy release and fluctuations in temporal statistics tempts one 
to interpret fracture via concepts of criticality and phase transitions. However, simple 
models have failed to convey the phenomenology and the origins of observed scaling 
laws are lacking. In this work we have performed extensive set of tensile experiments 
on paper sheets in order to discuss the concept of criticality in mode I loading with 
imposed strain. [TJ [2] . 

Typical stress-strain curve in tensile fracture of paper samples exhibits some non- 
linearity before the maximum stress, a c . This is a result of both plastic, irreversible 
deformation and to a much lesser degree of loss of elastic stiffness, or damage 
accumulation [3]. We study the damage accumulation by using "acoustic emission" 
(AE). This is the release of elastic energy due to microcracking on various scales ranging 
from far below the fiber size to, perhaps, millimeters as in the individual advances made 
by a notch in a tensile test. AE studies have been done in a wide variety for materials 
in science and engineering [4]. The common features in most are that the material 
failure process is complicated and disorder is present in the structure. Numerous studies 
have elucidated the statistical laws that describe AE. In general, these relate to usual 
mode I or mode III -type loading conditions, and the common feature is that they 
exhibit scalefree features. The probability distribution function (pdf) of event energies 
follows most often a power-law, with exponents in the range (5 = 1 ... 2. The event 
intervals are, similarly, found to obey such fat-tailed pdf 's. [3 El [TJ |2j [8] 

Here we consider paper as a test material or case of critical divergences during 
the fracture in the presence of structural randomness. This idea of criticality can 
be formulated in various ways, but essentially it implies the presence of a finite time 
singularity, so that a quantity such as the AE energy release during the approach of the 
failure becomes a function of l/(t c — t) where t c is the lifetime. In the following, we focus 
on the non-stationary AE signal from tensile experiments with imposed strain and study 
energy and event statistics of AE near the failure. We study in detail local averages 
of e.g released acoustic emission energy when t c , or critical point, is approached. The 
relation between critical point t c and the time at the maximum stress a c is discussed. 

The total number of observed acoustic emission events in all experiments is 800921. 
The high number of events in a single experiment (few thousands of events, comparable 
to rock and concrete fracture experiment cases [HI [10]) combined with the large number 
of samples allows us to study the non-stationarity of the acoustic emission signal. 

What kind of AE signature is produced by driving the fracture, e.g. by imposing 
a constant strain or creep stress, is not well understood. Thus, understanding the 
statistical properties of the acoustic emision and fracture precursors could provide a 
way to distinguish between different modes of loading. One possible application of this 
is in creating analogies with experimental findings in seismicity [TT] . 

The main result of the analysis is that we can distinguish two fundamental cases: 
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the dynamics in samples without a dominating notch and dynamics in the presence of a 
large defect with slow, stable crack growth. In the former case, there is evidence of such 
a divergence, upon approach of a sample- dependent critical time. This takes place after 
the sample maximum stress is reached. A qualitative explanation is perhaps offered by 
the dynamics which stems from the nucleation and growth of a dominating microcrack. 

To discuss the divergence, we consider simple models of statistical fracture. We 
show that variants of a minimal quasi-static fracture (random fuse network, RFN) 
model do not exhibit such features as the data. It is well-known that the RFN and 
other models do not reproduce quantitatively e.g. the exponents that characterize AE 
statistical properties. It also turns out that the divergence is not found in such models, 
however a RFN variant with gradual failure gets slightly closer [12] . 

The structure of the rest of this paper is as follows. In Section 2, we discuss the 
background and details of the experiments. Section 3 presents various aspects of the 
data analysis, and is finished with a subsection on modeling attempts. Section 4 presents 
discussion and summary. 

2. Experiments 

Normal copy paper samples were tested in the cross direction on a mode I laboratory 
testing machine of type MTS 400/M. The deformation rate e was 10 mm min~l. The 
AE system consisted of a piezoelectric receiver, a rectifying amplifier and continuous 
data-acquisition. The data-acquisition was free of deadtime. The stress was measured 
simultaneously directly from MTS 400/M using same time- resolution as AE-signal; thus 
the accuracy of time synchronization of the acoustic emission and stress-strain curves 
was below 10 /xs. 

During the experiment, we acquired bi-polar acoustic amplitudes simultaneously 
on two channels by piezocrystal sensors, as a function of time. Two transducers were 
attached directly to paper without a coupling agent. Each channel has 12-bit resolution 
and a sampling rate of 312 000 Hz. The transmission time from event origin to sensors 
is in the order of 5 /is. The acoustic channels were first amplified and after that held 
using sample-and-hold circuit. 

70 samples had initial notches of size 10 mm to achieve stable crack growth and 
100 samples were intact. Sample geometry was 150 x 50 mm 2 . Samples were tested in 
standard conditions: temperature of 22 degree in Celcius and relative humidity of 30%. 

The acoustic time-series was post-processed after the measurement by detection of 
continuous and coherent events, and the calculation of an event energy E^ is done as 
the integral of squared amplitudes within the event: E% = f A 2 (t')dt'. The event arrival 
time ti was taken from an instant when the amplitude raises above the threshold level. 
Dynamic range of the measurement device was 54 dB. The ensuing discrete set of events 
is characterized by set of pairs: {(t\, Ex), (i 2 , E 2 ), . . .}. Simultaneously, we measured the 
stress signal a(t) and, on the basis of that, we can define quantities also as a function 
of stress er(t). 
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Figure [T] shows the measurement setup and a single AE event. Time vs. stress 
curves from intact samples are shown in figure [21 which presents sample-to-sample 
variation in different experiments. The curves are shown both after scaling the sample 
strength (max. stress) and corresponding time to unity and in the usual way (inset). 
Both modes illustrate the sample-to-sample variation, and also in the proximity of the 
maximum stress. Figure [3] depicts time-stress-curves with a notch. As it is obvious, a{t) 
is of different shape and indicates stable crack growth. 

3. Results 

3.1. Temporal characteristics of the AE 

We start our analysis from a time ordered set of time-event pairs from an experiment: 
{(ti, Ei), (t 2 , E 2 ), (t n , E n )}. We divide the time interval (ti,t n ) into windows of length 
At and compute the event energy rate E(t) and event count rate n(t) for acoustic 
emission in each At window defined according to equations [I] and [2j 



where H(x) is the Heaviside step function. 

After energy and event count rates have been obtained, these quantities are averaged 
over all experiments. For the averaging, we have to use a similar choice for the time 
axis (or stress) and there is no natural choice. The time t = corresponds to triggering 
from the beginning of the tensile experiment, which is unphysical for any averaging. 
Before computing the average, we define origin of time and shift the signals n(t) and 
E(t) accordingly. 

As a physical beginning of the experiment, the origin of time is defined as where 
we find the maximum slope of time-stress curve (see figure [2]) and the intersection of 
the time axis and a straight line fitted to the point of maximum slope. This choice 
removes non-idealities at the beginning of the time-strain curve which often arise to 
experimental limitations, e.g. small slack in the sample when it is inserted to clamps. 
Figure H] shows E(t), the integral of event energy E(t), averaged so that the time is 
shifted to the beginning of experiments. The notched case shows stable crack growth. 
The unnotched energy integral depicts the onset of the acoustic emission when the plastic 
regime starts, but since time-to-failure scatters from sample to sample, the meaning of 
the average becomes subtle. The energy integral in the unnotched case is similar to the 
result presented in [6] . The differences between the result in [6] and figure H] are due 
to the fact that latter does not include events after the stress has started to drop after 
maximum has reached, the larger scatter in failure times and the samples being more 
ductile. It is interesting to note that the data for the notched samples shows over two 
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orders of magnitude in integrated energy a power-law scaling E cum oc t 16 . Assuming 
stable crack growth with a constant FPZ size, this feature could be interpreted as 
the crack length time dependence. Note that there is roughly one order of magnitude 
of difference between the samples without and with notches in total energy release 
magnitude. 

Next we test for apparent criticality, which would imply the divergence of the event 
energy rate near a critical point that is, the finite lifetime of the sample (corresponding 
to maximum stress, or crack instability). Discussing this is meaningful only in the 
unnotched case, where one can see a rapid increase in the event energy near the 
maximum of the time-stress curve. Conventionally, the lifetime t c is defined as the 
time at the maximum of stress o~(t), that is t c |max(cx(t)). However we discover that the 
time at maximum stress t c is different from the time at event the maximum of energy 
and event count rate t' c \ma,x(n(t)) . Schematically these differences are shown in the Fig. 
[5j The sample to sample distribution of the quantity (t' c — t* c ) is seen in figure [6] which 
shows that there is clear difference time at maximum event rate being larger. In other 
words, t' c is found without exception after the maximum stress. 

When we average quantities over samples by shifting the critical time to the origin 
before the averaging, as shown in figure [7J we find that event count rate n(t c — t) appears 
to follow a power law according to Eq. ([3]), with the lifetime as £ c |max(n(t)): 

n(t c -t) ~ (t c -t)- A (3) 

where we find the value A = 1.4 ± 0.1. The error is defined using the min-max method. 

Equation (j3j does not apply for the event energy release rate E(t c — t). A power 
law behaviour for the event count rate but not for the event energy implies a change 
in the average event energy during the experiment. This is shown in the figure [TO] 
average event energy as a function of normalized stress cr/a max during the experiment. 
The result shows that there is rapid increase in the average event energy near a c . We 
can also look at event and energy rates as functions of normalized stress (cr c — cr)/o~ c as 
shown in figure El The cumulative event count shows over a range of rescaled stresses 
an apparent power-law dependence. The cumulative event energy behaves in a way 
similar to what was found by Garcimartin et al [15] . but a power-law dependence is 
not clear, when we have the cumulative count as a reference.. The comparison is not 
straightforward due to different loading modes. The connection is further explored in 
the Discussion. As expected (on the basis of equation (j5J) the event rate does not seem 
to scale here in any particular fashion. 

3.2. Probability distributions of AE 

Next we study statistical distributions of event energies Ei and inter-event times. Inter- 
event times in an experiment are defined as set of Tj = U + i — tj. These exhibit almost 
without exception power-law statistics such that one defines probability distributions 
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P(t) and P(E) for the inter-event /waiting times and energies which are characterized 
by power-law exponents a and /3, respectively according to equations H] and [5j 



We find values (3 = 1.4 ± 0.1 and a = 1.3 ± 0.2 for event energies and inter-event 
times respectively. However, interpretation of these statistical distributions is difficult 
since a typical tensile test (or fatigue or creep test) is not stationary. The internal state 
of the material changes along the whole test duration, and statistical quantities such as 
P(E) or P(r) integrate over the whole history of a sample. Thus we study how these 
statistical distributions change e.g. in relation to t c or event count rate n. 

In figure [9] we show waiting time distributions from tensile experiments as a function 
of the event rate n. The event rate n is computed in a window. Windows are divided 
into different classes on the basis of an event rate. The data sets in the figure represent 
waiting time distribution in the event rate class. The data set label indicates the 
averaged event rate in units of s" 1 in the event rate class. Event rate classes are 
defined so that each class contains approximately 20 000 events. The difference between 
panels is the algorithm which identifies events: algorithm used in the lower figure allows 
overlapping of AE events while upper does not. The implication is that the algorithmic 
details are not important and differences we see in distributions is not due to the 
clustering of events at larger event rates. From waiting time distributions, it appears 
that there is a complicated change of the waiting time distribution P(t) as a function 
of event rate though one must be careful as the window size is still (in the time-domain) 
relatively large compared to t c . Perhaps, one can see signs of two different power- laws. 
It is interesting to note that the purest power-law behavior is obtained at small event 
rates (small strains) while for larger ones it changes to a more complicated distribution. 

The energy distributions are shown in figure [TTJ Due to large variation in the 
average event size (see figure [TO]) near time at maximum stress, we studied changes in 
the probability density P(E) close to t c . The distribution in the Fig. [TTJ contains events 
up to the time t c — t'. The time shift t' is shown in the label of the data set. The 
number of events is 6448 if we choose all events up to 10 seconds before t c , and 162425 
for all cumulated events up to 1 second after t c . Despite the large increase in the average 
event energy, overall form of the probability density function does not change as much 
as P(t) chang function of the event rate. Large event energies are cut off, when 

one considers the distribution which contains events up to the time which cuts off events 
after the event size has started to increase. 

3.3. Theoretical models 

In this chapter we compare experimental results to a random fuse network with a residual 
conductivity (RRFN). The model is analogous to the one proposed by Duxbury and Li 
[T2] . It is a minimal quasi-static model for fracture, which includes stress enhancements 



P{E) ~ E-f* 



(4) 
(5) 
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and residual bonding, leading to apparent plasticity and crack arrest. The model is an 
electrical analogy of an elastic lattice. The main motivation for using RRFN is that 
usual RFN simulations are much further away from the experiments. 

In the model we connect a set of fuses to a square lattice with a size L x L. 
Fuses are labeled as index j and associated with a constant conductivity gj = 1.0 
and with randomly distributed critical currents i C j. The distribution is uniform and 
characterized by its width: i c - G [1 — W, 1 + W]. For the breaking process we apply a 
slowly increasing voltage across the lattice until the hottest fuse approaches its critical 
current i C j. The important ingredient over the usual RFN is that fuses break at two 
phases: at the first phase its conductivity drops from gj = 1.0 to a residual conductivity 
gj = r. The fuse remains non-broken in the network with a residual conductivity until 
its current threshold is exceeded again and then its conductivity drops to zero, which 
finally corresponds to breaking the fuse and generating an event. We apply this rule to 
the conductivity of the hottest fuse and start increasing voltage from zero again. The 
process of applying conductivity drops and re-evaluating the hottest fuse is repeated 
until the conductivity of the network goes to zero and burned fuses form a percolating 
path across the lattice. 

When the residual conductance approaches unity, the model is identical to the 
perfectly brittle random fuse network. When the residual conductance is close to zero, 
stress enhancements create very tough fuses which are capable of arresting cracks and 
leading to plasticity. Width of the disorder parameter W leads to competing effect of 
the disorder and stress enhancements. 

By finding Vl-characteristic of the whole network, we obtain an entity which 
corresponds experimental stress-strain curves. We record a voltage-current pair (V(, I!) 
when a fuses conductivity drops to zero. We apply voltage-control for the ordered 
set {(V[, I'j), ...} by requiring that when the voltage V( overcomes the highest value 
previously recorded we get one point for the Vl-curve (Vfc, where k identifies an 
"acoustic emission" event. In the Fig. [12] we show all voltage-current pairs V(, V( from 
one numerical experiment and the envelope corresponds to the Vl-curve {(Vfc, Ik), ■■■}■ 
In the Fig. [13] are voltage-controlled Vl-curves from 99 simulations for the system size 



The event count rate is defined similarly as in experiments. We divide voltage-axis 
AV sized windows and compute number of events in a window according to Eq. [6] V max 
is the voltage when the current is in its maximum. The "acoustic emission" event in this 
model is associated to second conductivity drop of the fuse, i.e. when the conductivity 
drops to zero. Creating an event from all conductivity drops was also calculated, but 
this did not lead to an acceleration of the event rate. 



We performed set of simulations from system sizes 30 x 30 up to 124 x 124. The 
combination of disorder parameter W ~ 0.8 and residual conductivity r m 0.2 led 
to similar Vl-characteristic as in experiments: activity started after the experiment 



124. 




(6) 
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approached half of the critical voltage and we observed activity after the maximum 
voltage was approached (Fig. [T3|) . This implies that even large cracks are arrested 
before a dominating crack is formed. 

Fig. [TH shows the event count rate as a function of V/V max averaged over 99 
experiments. The event count rate starts to increase after 50% of the maximum voltage 
is approached. The event rate is in its maximum near V max . There is a small tail in the 
event count rate after the maximum stress, which becomes steeper when the system size 
increases. The increase in the event count rate is neither a power law nor critical. In 
the Fig. [15] there is the histogram of the maximum event rate a a function of V/V max . 
This depicts the difference between experiments and simulations: the maximum event 
rate occurs before or at the maximum current in simulations while in experiments it is 
observed after the maximum stress. However, there are still a few occurrences of the 
maximum event rate after the maximum current in simulations. We conclude that the 
time at the maximum event rate and the time at the maximum current are not identical. 
However, the relation is inverse to the experiment. 

4. Discussion 

We have compared fracture precursors using different "critical points" t c |max(<7 c ) and 
t' c |max(n) in mode I loading and imposing a constant strain rate for copy paper samples. 

A divergence is observed in samples without a notch, using time as a control 
parameter and looking at the event rate when approaching time at the maximum event 
rate. We have shown that the characteristic behaviour of the event rate E is different 
from the event rate h. The behaviour of the event rate h(t c — t) ~ t~ A might be taken 
to imply criticality, when t c is chosen to be the time at the maximum of the event rate. 
We note that the maximum of the event energy and event rate corresponds to a time 
which is observed after the maximum of the stress o c . 

A plausible argument to explain the causal nature of the AE observations (without 
a notch) is as follows. At the maximum stress it follows that d e a = E e e + E becomes 
negative. Thus E € is negative, which could be taken to indicate that one of the 
microcracks dominates and has started to grow. The event rate increases with the 
crack growth velocity, till the crack instability takes place. The divergence of the event 
rate could be then arise from a divergence of the crack velocity. The relation of the 
rate divergence and the concomitant increase of the event energy release presents a 
complicated problem. The rate of energy release rate should (recall we assume here a 
single propagating, stable crack) to be proportional to the new fracture process zone 
(FPZ) area created (if the AE event energy is related to the new FPZ area multiplied 
by fracture toughness). However, we have shown that the diverging quantity seems to 
be the event rate and not the event energy, in such strain imposed loading. 

Earlier experimental studies have shown indications of an event energy divergence 
when one imposes a constant pressure rate to a heterogenous material [T4"j IT3| IT5] . The 
critical exponent 7 in E ~ (e^p?} 7 was found to be 7 = 1.4 in [TJ], 7 = 1.27 in [13] 
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and 7 = 1.0 in [15]. In the case of an imposed constant strain rate there was not found 
to be any critical divergence of the energy release rate |16j . The analysis however was 
not done in the same spirit, and excluded the mechanism suggested above. Finally it 
was suggested that the real control parameter is time, since both imposing constant or 
cyclic stress seems to imply a critical divergence of the event energy [T7]. Our results 
include a somewhat similar behavior for the event energies as the maximum stress is 
approached. 

Statistical distributions of event energies P(E) and inter-event times P(r) have 
been measured earlier several times from paper. Power laws have been found with 
exponents /3 = 1.4 ± 0.1 and a — 1.3 ± 0.2 for event energies and inter-event times, 
respectively. We may contrast the current result to earlier work by Salminen et al. 
where acoustic emission has been measured from paper in the machine direction [6]. 
The energy exponent (3 is here slightly larger, compared to (3 = 1.25 ±0.10, which might 
be attributed to the much more ductile nature of paper as a material when stressed in 
the cross-machine direction. The waiting time exponent a is almost identical when it is 
integrated over whole experiment. 

In any case, in general most of the statistical signals that have been explored 
experimentally are still awaiting for theoretical explanations. For e.g. P{E), brittle 
fracture models allow to derive power-law -like scaling forms but with exponents ft that 
are in general far too large. To see if qualitative agreement could be found we compared 
the event rate from experimental data to the residual random fuse network. The event 
rate was found to increase when the voltage at maximum current was approached - 
contrary to perfectly brittle ordinary RFN. The result was not power-law divergence as 
in the experimental data. Note that introducing fracture toughness (residual fuses) to 
the model caused a change to the behaviour of the form of the event count rate: rapid 
increase of the event energy at V max was not found in the perfectly brittle RFN model. 
Further studies of the RRFN would seem to be of interest. Note that Amitrano and 
Helmstetter [18] have proposed a numerical model in order to model time-dependent 
damage and deformation of rocks under creep. In the model the 2D finite element 
method is used and separate time-dependent and time-independent damage progression 
and time-to-failure laws are introduced. Time-independent damage progression on an 
element is introduced as a gradual drop of an elastic modulus if the stress threshold was 
exceeded: Ei(n+1) = Ei(n)(l—D ) where Ei(n) is elastic modulus of an element after n 
damage events and D is a constant damage parameter. Quasistatic stress redistribution 
may induce an avalanche of damage events during a single loading step. The damage 
progression law and quasistatic stress redistribution is similar to conductivity drop and 
re-evaluation of the voltage in RRFN. The difference is in the damage law of an element. 
The result of the model is a power-law acceleration of the strain rate e and the rate of 
damage events h near the failure time t c . However, if we look at number of events 
as a function of strain j ~ ^ _^ 7< , the rate of damaged events decreases when t c is 
approached. The result is qualitatively in agreement what is observed on the residual 
random fuse network when it is assumed that every conductivity drop creates an acoustic 
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emission event. The gradual failure does not explain the event rate acceleration near 
the failure in the experiments. 

Eq. ([3]) is in essence identical to the Omori's law of aftershock rates. In the 
stationary paper peeling experiment one observes Omori's law with exponent close to 
A = 1.4 and it is also comparable e.g. to tectonic seismicity. In that case the generic 
temporal properties are further discussed in Ref. [IT]. Approaching a "main shock" in 
stationary experiment might be identical to approaching time at the maximum event 
rate in strain-driven non-stationary fracture experiment when we consider the behaviour 
of the event rate n. Fundamental differences of non-stationary and stationary AE signals 
are the different forms of the inter-event time distribution and the increase in the average 
event energy before failure (or main shock). However, this acceleration of the event 
rate in the stationary experiment is significant when averaging over a large number of 
sequences and the similarity may be coincidental. 

The limitation of the experimental setup is the limited choice of loading modes 
and lack of feedback control of the imposed strain or the stress. For example, driving 
material based on the acoustic emission event rate, as applied for the rock fracture in 
Ref. |19| . we could obtain an experimental test to the criticality and diverging quantities 
of the paper fracture; it would be interesting to see if it is possible to alter the event rate 
divergence near t c by controlling the strain based on AE feedback during an experiment. 
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Time t [s] 

Figure 1. The measurement setup: sample attached to MTS tensile testing machine. 
In this Figure single piezocrystal sensor is attached to back of the sample (actual setup 
consisted two piezocrystals). The figure on the right shows single AE event signal [6]. 
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Figure 2. Time vs. stress curves from all samples without notch. The maximum 
stress is scaled to unity using a ja c . This figure depicts typical stress response and its 
variations as a function of time. The inset contains same data without scaling of the 
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Figure 3. Time vs. stress curves from all samples with a notch. The maximum stress 
is scaled to unity using a/a c . This figure depicts typical stress response of notched 
samples. Average over all time-stress curves is shown. 




Figure 4. Averaged cumulative event energies as a function of time from the beginning 
of the experiment. The data is averaged over 100 and 70 samples. The notched case 
presents a power law increase of acoustic emission energy. The unnotched highlights 
the onset of acoustic emission, when the plasticity starts. 
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Figure 5. Schematic differences of t c |max(er(t)) and t' c \max(h(t)) in a single 
experimental time-stress curve. Inset shows time-stress-curve and main panel is 
magnified to the end of the experiment. 




Figure 6. The distribution of t' c \max(n(t)) — t c \m&x(cr(t)). The result indicates that 
time at the maximum of the event rate is larger than the time at the stress maximum. 
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Figure 7. The event energy release rate E and the event rate n as a function of 
t c — t. Energy release and event rates are computed in a window Ai=0.2s which is 
slidcd over AE time series. Critical time t c is defined as a time where the energy relase 
rate h reaches its maximum value. Event energy and event rates are averaged over 1 00 
experiments using AE from unnotched samples. The event rates decay as a power law, 
which is not true for event energies. 
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Figure 8. Cumulative energy E/Etot and cumulative event count n/ntot as a function 
of erc ~ a ■ Event energy and event counts are averaged over 100 experiments using AE 
from unnotched samples. 
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Figure 9. Waiting time distributions from tensile experiments. AE time series are 
divided to Ai=0.2 second time windows and waiting times r, and the event rate h is 
computed in a window. Windows are divided to different classes based on an event 
rate in the window and the data set in the figure indicates waiting time distribution 
in the event rate class. Label is the averaged event rate in the event rate class. 
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Figure 10. Average event energy, as a function of a ja c . This result shows that 
the event energy is not constant during the experiment. 
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Figure 11. The figure depicts evolution in the probability density of event energies 
when the time at the maximum stress t c is approached. The distribution contains 
events which are cumulated up to the time t c — t' . The time shift t' is shown in the 
label. The distribution becomes broader when the maximum event rate is approached. 
The solid straight line is a power-law fit P{E) ~ E~ 1A . 
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Figure 12. Voltage-current pairs from single RRFN simulation. System size is 
L = 124, disorder parameter is W = 0.8 and residual conductivity r — 0.2. Filled 
squares indicate Vl-pairs when a fuse is burnt, that is its conductivity drops to 
zero. The straight line corresponds linearly elastic behaviour. The envelope curve 
corresponds stress-strain curve of the numerical experiment. 
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Figure 13. 99 voltage controlled Vl-curves from RRFN simulations with system size 
L = 124. The disorder parameter is W = 0.8 and residual conductivity r — 0.2. The 
behaviour correspond qualitatively experimental stress-strain curves. 
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Figure 14. Event count rate near critical voltage from RRFN simulations. The data 
is averaged over 99 samples. Disorder parameter is W = 0.8 and residual conductivity 
r = 0.2. Increase in the event count rate is neither exponential nor a power-law. The 
tail of the event rate after the maximum current becomes steeper when the system size 
increases. The inset is a magnifigation of data with ten times smaller bin than in the 
main figure. 
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Figure 15. Histogram of maximum event rate occurrence as a function of V/V max . 
The result shows sample to sample variation of the relative location of the event rate 
maximum. Histogram is from L — 124 and it contains 99 samples. 



